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I. INTRODUCTION 


The preliminary definition phase included the following actions: 

1. Producing a Science Requiring Document-Draft 

2. Producing a Science Requirement Document- Preliminary 

3. Updating the Flight Program Proposal 

4. Project Review at NASA-Marshall Space Flight Center 

5. Research work as defined in 4.0 schedule, page 9 of the 
Statement of Work (May 25,1993) (see Fig 1.1) 

The first three items of this action plan have been delivered by the 
University of Alabama to NASA according to schedule. A Project Review Meeting 
was held at MSFC on June 29, 1993. Consequently, this part of the report will 
address the results of the research work performed in the Solidification 
Laboratory at the University of Alabama during the first six month of the project. 

The presentation of the experimental and theoretical work will correspond 
directly to the schedule given in Fig. 1.1. Tasks referred to in this report are El, 

E6, E8 and T1.1. 
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Figure 1.1 Schedule of the NASA/PEPSI project from Jan. '93 till Feb. '96. 




2. TASK El 

The scope of this task was testing and acceptance for operation of the 
Bridgman directional solidification furnace. 

The Bridgman directional solidification furnace was refurbished. This 
furnace can heat up to 1600 °C. Directional solidification can be performed in 
both the vertical and horizontal directions using this furnace. Some modifications 
were made to the initial design to increase its service life. Tungsten-Rhenium 
heating elements were replaced by Platinum-Rhodium elements to increase 
stability. The thermocouple connectors and power controllers were modified 
accordingly. Water connection to a pressure switch was introduced as a safety 
measure. The furnace is now ready for final testing. After the final testing, it will 
be interfaced with a computer to program the furnace operation. 


3. TASK E6 

The objective of this task was to design and build a resistance furnace 
with quenching capabilities that will be used to evaluate the local growth velocity 
of equiaxed grains. 

The design of the resistance furnace has been completed. This furnace 
uses molybdenum heating elements and it can heat to 1600 °C. It can work in 
vacuum as well as argon atmosphere. Quenching will be done by the drop- 
bottom technique. The sketch of this resistance furnace is given in Fig. 3.1. 

The furnace will be built after receiving all the parts, which are being 
ordered presently. 
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Figure 3.1 Sketch of the resistance furnace. 



4. TASK E8 

The objective of this task is to validate the model developed under task T 
2.1 using transparent organic materials. It must be noted that, since the physics 
of particle pushing by cellular interfaces is not clearly understood, it was 
necessary to perform the experimental task before the development of the model. 


4.1 Background 

Particle pushing or engulfment by advancing liquid/solid interfaces is 
relevant to a diversity of problems including particle chromatography, frost 
heaving, growth of semiconductor crystals, inclusions in steel, and metal matrix 
composites. The subject has received attention from experimentalists! 4 - 1 ' 4 - 6 ! and 
theoreticians alike.! 4 - 1 .4.3, 4.6,47-4.1 1] 

It is generally accepted that insoluble particles can either be pushed 
ahead by the interface, or engulfed in the solid phase. The behavior of the 
particles at the solid/liquid interface depends on the interface velocity. It has been 
demonstrated experimentally that there exists a critical growth velocity, below 
which the particles will be pushed by the interface, and above which they will be 
engulfed. f 4 - 1-4 - 7 ! Engulfment will normally lead to uniform distribution of particles, 
while pushing will result in particle segregation. If the interfaces converge on one 
another, as is the case in equiaxed solidification, the particle will be eventually 
entrapped in the solid at the end of local solidification. It is thus of great interest 
to characterize the critical velocity as a function of materials parameters and 
processing variables, to ultimately control particle distribution in the matrix. 

In a recently published paper, ! 4 - 12 ! an analytical model was proposed. It 
involved the following analysis: 

1) solving the temperature field for the particle/matrix configuration; 

2) calculating the shape of solid/liquid interface in the vicinity of the 

particle; 

3) calculating the forces acting on the particle; 

4) establishing the critical condition for the pushing/engulfment transition. 

The following equation was obtained for the equilibrium solidification 
velocity for the pushing/engulfment transition: 
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where d is the particle-interface distance, rj is the liquid viscosity, \i = Kp/Kl is the 
thermal conductivity ratio, Ds 0 is the interface surface energy, a 0 is the 
interatomic distance, Ap is the density difference between the particle and the 
molten matrix, r is the radius of the particle and n can take a range of values from 
1-7 depending on the various interactions (Coulombic, Van der Walls, etc.) 
between the particle atoms and the melt interface atoms. If gravity effects are 
neglected, which is reasonable for horizontal solidification, the critical velocity 
(maximum equilibrium velocity) is: 


V - a o Act 0 ( n - If 1 
c 3 r| p(n - 1) r l n / 


(4.2) 


Computer fitting of numerous experimental results I 4 - 3 ) produced an 
equation similar to Eq (4.3): 


V c = 



(4.3) 


where C is a constant and the exponent m ranges from 0.28 to 0.9. 

Equation (4.2) has been derived for the case of planar liquid/solid 
interface. However, under most practical growth conditions, the microstructure of 
alloys of commercial interest will be either cellular, dendritic or equiaxed. It is 
expected that the interaction between particles and non planar interfaces will be 
significantly different. 

Although some experimental studies have been carried out in metallic 
system with non planar interfaces I 4 - 13 ), a clear understanding of the 
particle/interface interaction is not yet available. The main difficulty is the 
impossibility of phenomena visualization because of the opacity of the systems 
used. This problem can be avoided by using organic transparent analogs of the 
metallic systems such as succinonitrile. 

Recently, some experimental work was conducted to examine the effect of 
particles on the morphological changes and on the nature of the entrapment 
processes in a transparent organic system. t 4 - 14 ) It was shown that the presence 



of a particle alters the solute gradient at the interface segment facing the particle. 
This causes interface perturbation. However, the effect of interface morphology 
on the pushing/engulfment transition was not documented. In another recent 
paper I 4 - 12 ], experiments were conducted with succinonitrile and SiC particles. It 
was found out that, while the behavior described by equation (4.3) is valid for 
larger particles, for small particles where V c is rather large, smaller values are 
obtained for the critical velocity than predicted by the power function. 

From the available information it is apparent that a more detailed analysis 
of the correlation between critical velocity and interface morphology is required. 


4.2 Experimental Procedure 

The matrix materials used in this work were succinonitrile (SCN) as the 
major component, and distilled water as the minor component. The following 
compositions were used: pure SCN (99% purity), SCN + 0.5 wt% water, and SCN 
+ 1 wt% water. 

The particles used in this study were polystyrene spheres with radii 
ranging from 0.5 to 7.5 pm. Directional solidification experiments were conducted 
in an experimental setup as shown in Fig. 4.1. It consists of a light microscope 
and a temperature gradient stage (TGS) as described by Ahuja t 4 - 15 ! The 
specimens were prepared by stirring particles into the molten SCN matrix, and 
then pulling the liquid in between the two glass slides through capillarity. The 
prepared specimen was placed between a hot plate and a cold plate separated 
by a distance of 6.67 mm. The temperature gradient at the solid/liquid interface, 
was changed by varying the temperature of the heater. For the experiments 
described in this paper, the temperature gradient was held constant at either 108 
K/cm (soak temperature T s = 85 °C), or 127 K/cm (T s = 95 °C). The translation of 
the TGS was adjusted through a speed controller. The translation velocity was 
varied between 1.5 and 17 mm/s. Particle behavior in front of the moving 
solid/liquid interface was observed and photographed with the microscope at up 
to 400X magnification. The particle size was measured using a stage micrometer 
and eyepiece reticule. 



microscope 



temperature 

controller 


Fig. 4.1. Experimental apparatus for directional solidification of transparent matrix 

materials. 


4.3 Results and Discussion 

A summary of the experimental results including the influence of 
composition, particle radius and temperature gradient on the critical velocity for 
engulfment is given in Table 4.1 . 


Table 4.1 Summary of directional solidification experiments with the SCN-H 2 O 
systems. 


Particle radius, 
mm 

Critical velocity, mm/s 

SCN 

G|_ = 108 K/cm 

SCN+0.5% H 2 O 
G[_ = 108 K/cm 

scn+i%h 2 o 

Gl = 108 K/cm 

SCN+1%H20 
Gl= 127 K/cm 

0.65 

16.80 

10.10 

6.31 

8.42 

1.25 

12.63 

10.10 

6.31 

8.41 

1.90 

ND 

8.42 

ND 

7.21 

2.50 

6.74 

6.31 

4.21 

5.61 

3.15 

3.37 

5.60 

3.88 

5.56 

3.78 

2.66 

2.50 

ND 

4.21 

4.37 

ND 

ND 

1.68 

ND 

4.41 

ND 

ND 

ND 

2.53 


ND: not determined 
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A Pure SCN. 

The experimental conditions were varied such that pushing or engulfment 
of polystyrene particles in SCN were observed. Fig. 4.2 shows a sequence of 
photographs demonstrating engulfment of a large polystyrene sphere (r = 5.6 
pm) in SCN solidifying with a plane front. The solidification velocity was 1.48 
pm/s and the temperature gradient was of 108 K/cm. It can be noticed that the 
liquid/solid interface is distorted by the approaching particle, and a depression 
forms. This type of distortion was observed in all the experiments conducted in 
this work. However, according to the proponents of the "thermal conductivity 
criterion" model [4.8,4.i6] f when p = % <1, which is the case for the system 

polystyrene particles - SCN used in this study, a bump rather than a depression 
should form at the liquid/solid interface. Indeed, a qualitative analysis of the 
thermal field by Bolling and Cisse t 4 - 8 !, as well as the analytical model proposed 
by Stefanescu et al. I4.6,4.i2] f suggest the formation of such bumps on the 
interface. Nevertheless, these analyses ignore the latent heat evolved during 
solidification. This heat can be "trapped" between the particle and the interface, 
because of the lower thermal conductivity of the particle. Thus the interface will 
be depressed, and a depression rather than a bump should form. Indeed, a 
numerical model proposed by Sasikumar et al. t 4 - 10 !, which includes the entropy 
of fusion, describes the formation of a depression at the interface even for p < 1. 

Furthermore, it is reasonable to assume that, since the "pure" SCN 
contained about 1% impurities solute was rejected and piled-up in the liquid film 
between the particle and the solid/liquid interface. This solute accumulation will 
result in a temperature depression at the interface, and thus may favor formation 
of a depression. 

For smaller particles pushing was easily observed. This is illustrated in 
Fig. 4.3 for two small polystyrene spheres having radii of 1.9 and 3.0 pm, 
respectively. The particles are pushed by the planar SCN solid/liquid interface at 
the imposed velocity of 5.5 pm/s. As evident from Table 4.1, this velocity is 
smaller than the critical value. A depression of height smaller than the radius of 
the particles is observed. 
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Fig. 4.2. A spherical polystyrene particle (r = 5.6 pm) being engulfed by the 
planar SCN solid/liquid interface at V=1.48 pm/s and G|_=108 K/cm. (a)-(d) 

increasing time. 



Fig. 4.3. Two polystyrene spherical particles (r=1.9 pm and 3.0 pm) being pushed 
by the planar SCN solid/liquid interface at V=5.5 pm/s and G[_=108 K/cm. 


The dependence of the critical velocity on particle radius is plotted in Fig. 
4.4. The full squares represent the data for velocities resulting in planar interface. 
The empty squares correspond to velocities at which cellular interface was 
obtained. Apparently, the critical velocity decreases with increasing particle size 
as predicted by theoretical models. The solid lines on the figure are the least 
square best fit through the data points with Eq (4.3), assuming that m = 1, as 




suggested by Eq (4.2). The corresponding equations, for velocity expressed in 
pm/s, are: 


for planar interface only V c = (4.4) 

for all experimental data V c = (4.5) 



Particle radius, pm 

Fig. 4.4. Correlation between the critical velocity and particle radius in pure 

SCN (G L =108 K/cm). 


It can be seen that when all data are included, the V c - r curve is moved slightly 
to the left. This means that the break down of the interface resulted in lower 
critical velocity than expected for a planar interface. Nevertheless, pushing was 
still possible. At velocities above the critical velocity, particles were incorporated 
through engulfment. 

These results compare favorably with the experimental findings of Korber 
et al. I 4 - 17 ! where V c was found to vary proportionally with r 1 for latex particles in 
an ice-water system. 

Taking n = 2 as proposed in ref. [4.6,4.9,4.12], Eq.(4.2) can be used to 
calculate the difference in surface tension, Aa 0 . Using constants from the 

Appendix, one obtains Aa 0 = 2.53xl(T 4 N/m. When this value is used a good 
curve fit is obtained for the experimental data in Fig. 4.4. 




B. Succinonitrile-water systems. 


The influence of solute on the critical velocity was studied with polystyrene 
particles in SCN containing 0.5 or 1.0 wt% distilled water. In the two-component 
systems the interface became morphologically unstable at lower growth velocities 
as compared with pure SCN, because of solute rejection and constitutional 
undercooling. Cellular structures were readily observed. Figs. 4.5 through 4.7 
show particle/interface interaction in SCN-0.5% water as the interface velocity 
was increased from 2.1 to 8.42 pm/s. 

In Fig. 4.5 it is seen that two small polystyrene spheres (r=1.9 pm and 3.0 
pm) were pushed by the planar solid/liquid interface with a velocity of 2.1 pm/s. 
The planar interface began to break down as the solidification velocity was 
increased to 5.05 pm/s. At this point, particles were pushed by a non planar 
interface (Fig. 4.6). The particle with the radius of 3.0 pm (particle on the left) was 
first engulfed by a non planar interface when further increasing the solidification 
velocity to 8.42 pm/s, while the particle with the radius of 1.9 pm (particle on the 
right) was still pushed by the cellular interface, and was eventually entrapped at 
cells boundaries (Fig. 4.7). 



Fig. 4.5. Particle/interface interaction in SCN-0.5 wt% water. Two polystyrene 
spherical particles (r=1.9 pm and 3.0 pm) are pushed by the planar solid/liquid 
interface at V=2.1 pm/s and G|_=108 K/cm. 
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Fig. 4.6. Particle/interface interaction in SCN-0.5 wt% water. Interface velocity 
was increased at V = 5.05 pm/s. The planar interface begins to break down. 
Particles are still pushed by a non planar interface. 



Fig. 4.7. Particle/interface interaction in SCN-0.5 wt% water. Interface velocity 
was increased at V = 8.42 pm/s. The particle with the radius of 1 .9 pm is pushed 
by the cellular interface, then entrapped at cell boundaries. 


Similarly, Figs. 4.8 through 4.10 represent the various phases of the 
solidification process in SCN-1 .0 wt% water when increasing the velocity from 2.1 
to 6.31 pm/s. Compared with the SCN-0.5 wt% water system, in the SCN-1. 0 
wt% water system, the planar interface began to break down at a lower interface 
velocity of 3.6 pm/s. In general, as expected, the velocity at which interface 
instability occurs decreases as the water content increases (Fig. 4.1 1). 


On Figs 4.8 and 4.9 it can be seen that some particles pushed by the 
interface are moving transverse to the interface, until they reach an intercellular 
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boundary. Then, these particles continue to be pushed at the cell boundaries. As 
the interface velocity is further increased to 6.31 pm/s, the interface becomes 
cellular, with subsequent particle entrapment (Fig. 4.10). 


Fig. 4.8. Particle/interface interaction in SCN-1 wt% water. Particles with radii of 
1.9 to 3.1 pm are pushed by the planar solid/liquid interface at V = 2.1 pm/s and 

G|_ = 108 K/cm. 


Fig. 4.9. Particle/interface interaction in SCN-1 wt% water. Particles with radii of 
1.9 to 3.1 pm are still pushed by a solid/liquid interface that has begun to break 
down at V = 3.6 pm/s and G|_ = 108 K/cm. 
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Fig. 4.10. Particle/interface interaction in SCN-1 wt% water. Particles with radii of 
1.9 to 3.1 pm are entrapped at cell boundaries by a cellular solid/liquid interface 

at V = 6.31 pm/s and Gl = 108 K/cm. 



Fig. 4.1 1. Experimental interface velocity above which the liquid/solid interface 
becomes unstable as a function of water content (Gl= 1 08 K/cm). 


The experimental values for the critical velocity as a function of particle 
radius for SCN-water systems are plotted in Fig. 4.12. It can be seen that, as the 
interface becomes unstable, the Vc-r correlation strongly deviates from the power 
function observed for pure SCN (compare with Fig. 4.4). A moderate deviation is 
observed when particles are engulfed by a non-planar interface. A strong 
deviation occurs when particles are entrapped by a cellular interface. In general, 
as the particle size decreases, and the velocity increases, particle incorporation 
in the solid occurs through entrapment rather than through engulfment. 





All the experimental results for samples solidified under a gradient of 108 
K/cm are plotted in Fig. 4.13. The full points represent as before the data for the 
experiments were the solid/liquid interface remained planar. The curve 
represents the best fit through all the points obtained for planar interface, that is, 
including SCN, SCN + 0.5% water, and SCN + 1.0% water. The equation 
corresponding to this curve is: 

= (4.6) 


entrappment engulfment 



Particle radius, pm 

Fig. 4. 12. Experimental variation of the critical velocity with particles radii 
in SCN/water systems (G|_=108 K/cm). Full data points correspond to 

planar interface. 

It can be seen that there is a minor difference as compared with Eq (4.4) for pure 
SCN. This seems to indicate that the solutal field plays a minor role on the critical 
velocity, within the experimental conditions. Nevertheless, it may be inferred that 
the solute slightly decreases the critical velocity. Deviations from the law in Eq 
(4.6) are stronger as particle size and interface stability decrease. 

From the well-known constitutional supercooling criterion proposed by 
Chalmers et a/.I 4 - 18 ! one can obtain the velocity at which interface instability will 


occur: 
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where Gl is the temperature gradient in the liquid at the interface, C 0 is the alloy 
composition, k is equilibrium partition ratio, D|_ is diffusion coefficient of solute in 
the liquid, and mL is the slope of the liquidus line. Thus, as the temperature 
gradient increases, the velocity for instability increases. Indeed, as shown on Fig. 
4.14, as the temperature gradient was increased from 108 to 127 K/cm, the 
velocity for instability (represented by full horizontal lines) increased from 3.6 to 
4.6 pm/s. As a result of this change the deviation of the V c - r curve from the 
power function occurs at higher velocities. In addition it is seen that for the SCN + 
1 wt% water system an increase in the temperature gradient at the interface from 
108 to 127 K/cm results in a systematic increase in the critical velocity. This 
increase seems to be larger for the case of cellular interface (velocities above the 
limit of stability), than for the case of planar interface. The V c increase with 
temperature gradient for the planar interface is in agreement with some 
theoretical predictions. Indeed, the numerical model proposed by Potschke and 
Rogge I 4 - 9 ! for systems with impurities suggested the following dependency: 
V c ~G L °- 5 r _1 . The gradient dependency of V c is due to the fact that the 
interface curvature is greatly affected by the temperature distribution in the 
system. 



Particle radius, pm 



0 1 2 3 4 5 6 

Particle radius, pm 


Fig. 13 Correlation between the critical pjg_ 14 Experimental variation of the 

velocity and particle radius in pure critical velocity with the radius of the 

SCN and SCN + various amounts of particles for different temperature 
water (Gl=108 K/cm). gradient (SCN + 1 % water). 




4.4 Conclusions 

The role of interface morphology on the critical velocity for particle 
engulfment in SCN, and its mixtures with water (0.5 and 1 wt%) was investigated 
through directional solidification experiments. Polystyrene spheres with radii 
ranging from 0.5 to 7.5 pm were used. 

It was found that the experimental correlation between V c and particle 
radius for planar interface in pure SCN can be described through the relationship 
V c - - 1 - 3 l 5-. This agrees well with theoretical predictions. When all the 

experimental data for planar interface, including SCN, SCN + 0.5 % water, and 
SCN + 1 % water were used the correlation changed to V c = This suggests 
that the solute slightly decreases the critical velocity. When the data for cellular 
interfaces were also included in the analysis, the correlation became V c = . 

This means that the break down of the interface further decreased the critical 
velocity. In fact, as the interface became unstable in the SCN + water systems, 
the V c -r correlation strongly deviated from the V c = cons./r law, which was found 
to describe well all experiments with planar interface. 

It was also found that the critical velocity increased with the temperature 
gradient for both planar and cellular interfaces. No evidence of depression 
formation at the liquid/solid interface was observed, although the thermal 
conductivity of the particle was smaller than that of the matrix. This apparent 
discrepancy with the theory is attributed to the fact that some models ignored the 
role of the latent heat evolved at the interface. 


4.5 Appendix 1: Numerical Values Used in Calculations 

Viscosity of succinonitrile h = 2.59x10*3 Kg m*is*i I 4 - 19 l 
Atomic distance in succinonitrile, a Q = 5.9x10' 10 m 
Thermal conductivity: 

for succinonitrile: K|_ = 0.223 Wm^K' 1 1 4 - 19 ! 
for polystyrene: Kp = 0.0794 Wm^K* 1 1 4 - 19 ! 
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alloy composition 
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atomic distance 
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diffusion coefficient of solute in 
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temperature gradient in the 

dc 

critical distance between the 

liquid at the interface 

particle and the interface 

Kl 

thermal conductivity (matrix) 

9 

acceleration due to gravity 

K p 

thermal conductivity (particle) 

k 

equilibrium partition ratio 

T s 

soak temperature 

m L 

slope of liquids line 

V 

interface velocity 

r 

particle radius 

Vjnst 

velocity at which interface 

u 

thermal conductivity ratio 
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h 

liquid viscosity 

V c 

critical velocity 
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density difference 
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equilibrium velocity 
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interfacial energy difference 
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5. TASK T1.1 

The objective of Task 1 is to further improve the existing analytical model 
for pushing/engulfment of particles during solidification of metallic alloys with 
planar solid/liquid interfaces. The scope of sub task T 1.1 is to develop a model 
that includes an analysis of both the thermal and force fields. 

An enhanced analytical steady state model for predicting the critical 
velocity of engulfment of an insoluble particle by the advancing melt interface 
was developed. This model couples the thermal and force fields acting on an 
insoluble particle at the liquid -- solid interface. Equations that determine the 
equilibrium distance for the force and thermal fields as a function of velocity 
were derived. The solution of these equations gives a unique value for the 
critical velocity of engulfment, and for the critical particle -- interface distance 
during pushing. The role of the thermal gradient, of particle and matrix thermal 
conductivities, and of the interface force is discussed. The model is validated 
with experimental results for succinonitrile and SiC or polystyrene particles. 
Literature data on the critical velocity for aluminum based composites containing 
alumina or SiC particles, and for the water -- tungsten systems also compare 
favorably with calculated results. The model is applicable to particulate metal 
matrix composites. 


5.1 Background 

One of the processing methods for particulate reinforced metal matrix 
composites (MMC) consists of mechanical mixing of particles in the liquid alloy. 
In the solidification interval the structure of such a composite includes three 
phases: liquid alloy, solid particulate reinforcements, and solidified alloy grains. 
The as-cast microstructure is the result of solidification kinetics that governs the 
growth of the equiaxed grains, and of particle/grain interface interaction that 
governs pushing or engulfment of particles into the growing grains. As a result of 
these two simultaneously occurring phenomena, the reinforcing particulates can 
be distributed intergranular, intragranular or mixed at the end of solidification. 
Since this distribution determines the properties of the composite its prediction is 
of considerable interest. 



In a first approximation the issue of particle/grain interaction can be 
treated in a similar manner as the particle/planar interface interaction, assuming 
that the grain is large as compared with the particle. It is generally accepted that 
when a moving solidification front intercepts an insoluble particle, it can either 
push it or engulf it. In the case of engulfment the solid grows around the particle. 
It has been demonstrated experimentally I 5 - 1 ! that, at a given critical velocity, 
fixed for the system, a transition from particle pushing to particle engulfment 
occurs. If the solidification interface becomes unstable resulting in formation of 
dendrites or equiaxed grains, two or more solidification fronts can converge on 
the particle. Even if the particle is not engulfed by one of the fronts, it will be 
eventually entrapped in the solid at the end of the local solidification. The 
problem is of rather general interest to such problems as particle distribution in 
ceramic particulate reinforced metal matrix composites, growth of monotectics, 
inclusion segregation in castings, frost heave in soils, and emulsion of organic 
cell suspensions in ice in cryobiology. More recently a novel theoretical approach 
to the mechanism of heterogeneous nucleation of pores in metallic systems has 
been developed on the basis of the particle pushing theory l 5 - 2 !. 

The interaction between a particle and an advancing planar liquid/solid 
interface has been discussed by a number of authors 15 . 3 - 5 . 6 ] An analytical 
model has been proposed by two of the authors of this paper l 5 - 7 ) (referred to as 
the SDKM model), and further refined recently (the SAS model) I 5 - 8 !. 

This last model involved the following analysis: (1) solving the temperature 
field for the particle/matrix configuration in order to calculate the shape of the 
solid/liquid interface in the vicinity of the particle; this is needed because the 
interface will locally deviate from planarity if the thermal conductivity of the 
particle, Kp, differs from that of the matrix, K|_, while a macroscopically planar 
interface may still be maintained; (2) solving the force field acting on the particle; 
(3) establishing the critical condition for the pushing/engulfment transition. The 
solidification latent heat was ignored in the solution of the thermal field. The 
model considers the behavior of a particle in front of a macroscopically planar 
solid/liquid interface; however, local deviations from planarity are allowed (Fig. 
5.1). 
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a) b) 

Fig. 5.1 Schematic drawing showing (a) forces acting on the particle and 
coordinate system for calculation of the force equilibrium distance, dF, and (b) the 
coordinate system for calculation of the thermal equilibrium distance, dj, for 

Kp/K m > 1 . 


The force field acting on the particle was analyzed based on Uhlmann et al. 
model f 5 - 3 !. This model assumes that the interface energy changes with the 
particle/interface distance as shown in Fig. 5.2. 



Fig. 5.2 Variation of interface energy with particle/interface distance l 5 - 3 ). 



If one writes: 


o 0 = a sp - (o s i + 0| P ) (5.1) 

then, the interface energy change as the particle approaches the interface can be 
written as: 


Ao = Ao 0 [ — ^-rl" (5.2) 

\a 0 + dj 

The subscripts s, I, and p stand for solid, liquid and particle, respectively. 
According to Uhlmann et al. I 5 - 3 ) n may have values of the order of 4 or 5. It must 
be noted though that, if Coulombic forces are considered, n = 2, while for Van der 
Waals forces n = 7. Coulombic forces are active over distances of the order of 
10' 8 m, while Van der Waals forces are effective at distances of the order of 10' 
10 m. Unfortunately, the actual nature of the interaction is not known. Overlapping 
interactions may result in any values of n, including n < 2. Accordingly, it will be 
assumed that values of n between 1 and 7 are conceivable. 

From the force balance on the particle (the drag force Fd and the 
interaction or repulsive force Fj n t), and neglecting gravity effects (which is a good 
approximation during horizontal directional solidification), the following equation 
was obtained for the equilibrium solidification (interface) velocity for the 
pushing/engulfment transition: 

v _ Aap d F I — ao — ( 5 .3) 

3 r| p R la 0 + dF/ 

where Ds 0 is the interface energy at minimum separation, dF is the steady state 
particle/interface distance for the force field, h is the liquid viscosity, m = Kp/K|_ is 
the thermal conductivity ratio, R is the radius of the particle, and a 0 is the 
interatomic distance. A detailed discussion of the nature of the forces acting on 
the particle, as well as the mathematical derivation of Eq (5.3) was provided in 
Ref. I 5 - 8 ). 

For n = 2, equation (5.3) gives two solutions for dF, for a given interface 
velocity, V e . In order to remove this uncertainty, in the SAS model it was 
assumed that the initial velocity corresponds to the unique value of dF, which is 
obtained through maximization of V e with respect to dF- The physics behind this 



maximization can be understood from Fig. 5.3, which shows the variation of V e 
versus d calculated for an Al-SiC system. As the interface velocity is increased 
the equilibrium distance at steady state decreases until a critical distance is 
reached for the maximum equilibrium velocity, which is considered to be the 
critical velocity. This is done since a further increase in velocity does not have 
any solution on the basis of the mathematics used in SAS-92. It is thus implied 
that at interface velocity higher than Vc, steady state conditions cease to exist, 
and the particle is engulfed. It must be noted that selecting the distance 
corresponding to the maximum equilibrium velocity to be the critical distance is 
only an assumption. To find a unique solution for Eq.(5.3) another equation is 
needed. The quest for an unique solution was thus the goal of the present 
research. 



0.0 0.1 0.2 0.3 0.4 

d, jim • 0.01 

Fig. 3 Calculated correlation between the equilibrium steady state interface 
velocity and particle - interface distance for an Al - SiC system 


5.2 Theoretical Work 

A. Analysis of the Force Field (SDKM /SAS Models) 

Eq.(5.3) was plotted in Fig. 5. 4 for an Al/SiC system for various values of 
n. Particle radius was considered to be 50 mm, and the other constants used for 
calculation are given in Table 5.1 in the Appendix of chapter 5. It can be seen 
that, with the exception of the case n = 1 , for each velocity there are two solutions 




for dF- Obviously, only those solutions giving dF > a 0 can be accepted, since dp < 
a 0 means particle engulfment. Thus, at steady state, for each interface velocity a 
minimum equilibrium distance between the particle and the interface will exist. In 
other words, as the velocity increases from a value smaller than the critical 
velocity, V cr , the particle -- interface distance decreases to maintain steady state. 
If no other effects are considered, this distance can decrease until the maximum 
of the V -- dF curve is reached. If the velocity increases above the maximum 
value for dF, steady state does no longer exist, and engulfment occurs. In the 
SDKM model I 5 - 7 ) it was arbitrarily assumed that dcr = a 0 . In the SAS model I 5 - 8 ) 
the maximum of the V -- dF curve was considered to be the critical distance, and 
the corresponding velocity -- the critical velocity. For n = 2, similar equations were 
thus derived for V cr for both models. 



Distance from interface, pm 

Fig. 5.4 Influence of assumption made for n in Eq.(5.1) on the critical distance 
obtained through maximization of the equation for Al/SiC composite. 


The locus of d C r as a function of n is shown in Fig. 5.4. It can be seen that 
values for n > 3 are unacceptable, since critical distances smaller than the atomic 
distance are predicted. 

In the SAS model mainly the force field was used to evaluate the critical 
velocity and the critical distance. The thermal field was only considered in the 
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calculation of the drag force, since it changes the shape of the interface. 
Nevertheless, its role in the establishment of the particle -- interface equilibrium 
distance was neglected. 

B. Analysis of the Thermal Field 

To find a unique solution to Eq.(5.3), that is to uniquely determine the 
critical velocity and the critical distance, an additional correlation between V and 
d is required. This correlation can be obtained if it is assumed that at steady state 
the thermal field will also impose an equilibrium distance for a given V. This is the 
thermal equilibrium distance, dy. As shown in Fig. 5.1, the interface will bend in 
the vicinity of the particle to accommodate the change in heat flux resulting from 
the Kp/K|_ difference. The change in the radius of curvature of the interface, and 
therefore in the particle -- interface distance, will continue until equilibrium is 
reached. 

The governing equation is the time independent form of the heat 
conduction equation with moving boundaries. A two- or three-dimensional 
solution will involve boundary conditions that are characteristic for the particular 
experimental setup used. Thus, only a one dimensional solution will be attempted 
at this point. The governing equations in Cartesian coordinates are: 

9^T 

for the particle: — = 0 (5.4) 

dz 2 

-) T 

for the matrix: ^ — = 0 (5.5) 

dz 2 a dz 

where T is the temperature, V is the solidification velocity, and a is the thermal 
diffusivity. It was assumed that the thermal conductivities of the solid and liquid 
are identical and are different from that of the particle. Heat transfer in the x and y 
directions was ignored. The temperature gradient in the melt, far away from the 
particle (see Fig. 5.1) is constant. The origin is in the center of the particle. The 
following boundary conditions were used: 


BC1: K S G S - K L G L = AH V 


(5.6) 
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BC2: At z = 0, T p = T soak - G L L where 

the temperature 

gradient in the liquid is defined as Gl = , Tsoa jL TiRt 

L + 2K + dj 

and since dj « R, can be approximated to:G L = Ts ? ak ~ p int 

L + 2R 

(5.7) 

«* &L-6L 

(5.8) 

BC4: (Tmjz = _R = (Tp)z = _r 

(5.9) 


where, Ks, K|_ are the thermal conductivities of the solid and liquid, respectively, 
G s is the temperature gradient in the solid, G L is the temperature gradient in the 
liquid film between the particle and the interface, Gl is the temperature gradient 
in the liquid, T soa k, T int , T m , and T p are the soak, interface, matrix, and particle 
temperatures, respectively, L is the length of the temperature field over which the 
gradient is measured, and AH is the latent heat of fusion. The soak or the set 
temperature is measured at a distance L from the solid/liquid interface. The 
solutions of Eqs. (5.4) and (5.5) are as follows: 

T m = A + B exp/-^-z) (5.10) 

\0tm / 

and T p = C + Dz (5.11) 

The constants A, B, C, and D are obtained by using the boundary conditions, 
Eqs. (5.6) through (5.9) in the solutions, Eqs. (5.10) and (5.11). They are: 

A = Tsoak - GlL -[Sm + B] (^M)exp[ (5.12) 

B Hv)( M tr Y ) exp (^ R+dT) ] <5 - 13) 

c = Tsoak • G l L - R 


(5.14) 



(5.15) 


D = 


( K S G S - LV | 
4Kl ) 


exp 


Vd T 
. OCm. 


Substituting the values of the constants A and B in Eq. (5.10) the final solution for 
the temperature of the matrix (liquid) at any position along the z axis is calculated 
to be: 


£«***■>] 


T m = T soak -G L L-^- 

M' * 



(5.16) 


To calculate the thermal equilibrium distance, it is now necessary to obtain a 
particular solution of Eq. (5.16) for a position on the solid/liquid interface. That is 
for: 


z = -(R+d T ) where T m = T in) 

Substituting in Eq.(5.16), the following solution is obtained for the thermal 
equilibrium particle - interface distance: 


dr^ln 


(T SO ak " T| n t 


Iv) 

|K s Gs - LVj| 

1 l 

(a m , R 
1 V p) 

(K S G S - LV 

l k l 

) I 


(5.17) 


Substituting the value of (T soa k - Tj n t) from Eq.(5.7) , Eq.(5.17) can be simplified 
to: 



(5.18) 


C. Coupled force and thermal fields 

It is assumed that the critical distance is that one at which the particle is in 
equilibrium in both the force and the thermal fields. The force field and the 
thermal field must thus be coupled, that is: 
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dF = dj (5.19) 

To couple the force field and the thermal field Eqs.(5.3), (5.18) and (5.19) must 
be solved simultaneously, graphically, analytically, or numerically. They are 
plotted in Fig. 5.5 for a succinonitrile/SiCp system for particles of radius 18 mm 
and different values of n. The data used for calculation are given in Table 5.1 
(Appendix chapter 5). a 0 was taken as the molecular length of succinonitrile. The 
intersections of the two equations give unique values for both the critical velocity 
and the critical distance. For example, for n = 1.5 a critical velocity of 4 mm/s and 
a critical distance of 8.74 mm are calculated for the 4 mm radius particle. It can 
be also noticed that while the equilibrium distance is practically unaffected by the 
value taken for n, on the contrary, the critical velocity is highly dependent on n. 
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Fig. 5.5 Influence of assumption made for n in Eq.(5.2) on the calculated critical 
velocity for the succinonitrile/SiCp system. Particle radius is 4 pm. 


5.3 Experimental Work 

To validate the proposed model a number of experiments were conducted 
on transparent organic material (succinonitrile) and various particles. The 




experimental setup for directional solidification of transparent matrix systems is 
shown schematically in Fig. 5.6. A complete description of the experimental 
apparatus is given in Ref. I 5 - 9 ). Horizontal solidification was preferred to eliminate 
the role of the gravity force on the particle -- interface interaction. 


microscope 



temperature 

controller 


Fig. 5.6 Temperature gradient stage for directional solidification of transparent 

matrix materials. 

Specimens were prepared by drawing a previously prepared mixture of 
particles and succinonitrile between two quartz slides. The slide was placed 
between the heater and the cooler. In this way a temperature gradient was 
imposed over the slide. The temperature gradient at the solid -- liquid interface 
was changed either by decreasing the distance between the heater and the 
cooler, or by varying the temperature of the heater. The solidification velocity was 
controlled through the translation velocity of the slide imposed by the motor. 

The matrix in all the experiments consisted of succinonitrile. Two kinds of 
particles were chosen for this study: polystyrene particles of diameters varying 
from 3 to 80 pm and SiC particles 8 and 36 pm in size. Polystyrene has a lower 
thermal conductivity than succinonitrile, while SiC has a higher conductivity. 

The interface velocity was set at 3.68 mm/s and the temperature gradient 
was measured in the vicinity of the melt interface. Typical measured data 
required in Eq.(5.17) are given in Table 5.1. 

In the experiments with SiC particles dispersed in succinonitrile (SCN) it 
was found out that at a solidification velocity of 9.2 mm/s, a large particle (R = 18 







mm) was engulfed (Fig. 5.7,a), while a small particle (R = 4 mm) was pushed. 
When the solidification velocity was increased to 21 mm/s, the small particle 
induced the formation of a trough on the interface, and was eventually engulfed 
as the borders of the trough closed behind it (Fig. 5.7, b). 


Table 5.1 Experimental data used in the calculation of the thermal equilibrium 

distance. 


Soak Temperature, T s (°C) 

64 

71 

75 

80 

Temp. Gradient, G|_ (°C/m) 

4.96*1 0 3 

6.24*1 0 3 

9.06*1 0 3 

1.08*10 4 


A similar experiment was run using 75 mm diameter polystyrene micro 
spheres dispersed in SCN. The thermal conductivity of polystyrene (0.0794 W m* 
1 K* 1 ) is lower than that of SCN. Under a solidification velocity of 2.5 mm/s the 
particle was engulfed. Again, trough formation occurred on the melt interface. 



a) large particle (small V cr ); V gr0 wth = 9-2 mm/s, time step 3 s, R = 18 mm 



b) small particle (large V C r); V gr0 wth increased from 9.2 to 21 mm/s, time step 4 s, 

R = 4 mm 


Fig. 5.7 Experimental observation of engulfment of SiC particles by the US 

interface in succinonitrile. 


5.4 Discussion 

Since the value of the Kp/Km ratio influences not only the shape of the 
interface but also the particle -- interface interaction, it is reasonable to divide the 
discussion in two parts, according to the value of the ratio. 

A. Case 1: Kp/Km > 1 

The physics behind Eqs. (5.3) and (5.18) and Fig. 5.5 may be understood 
as follows. Consider a particle in the vicinity of the interface moving at a velocity 
V > V cr (Fig. 5.8). The force field will allow the particle to approach the interface 
up to a distance dp < d cr and dF < dj (Fig. 5.9, a). Consequently, the thermal field 
will "push" the interface trying to reestablish the equilibrium distance dj (Fig. 
5.9, b). The distance becomes larger than the equilibrium distance for the force 
field. Consequently, the particle will again approach the interface to a distance dF 
(Fig. 5.9,c). This interplay between the two equilibrium distances will continue, 
and steady state will not be reestablished. The particle will sink continuously in 


the solid. The borders of the gap will close behind the particle, resulting in 
engulfment (Fig. 5.9, d). Coupled with this is the phenomenon of general 
destabilization of the interface because of the alteration of the thermal field. 



Distance from interface 


Fig. 5.8 Schematic representation of the difference between the equilibrium 
distance for the force field, dF, and that for the thermal field, dj, when interface 
velocity is higher than the critical velocity. 



Fig. 5.9 Trough formation and engulfment of a particle by an interface traveling at 
a velocity higher than the critical velocity. 




Experimental evidence of significant trough formation because of the 
influence of the particle on the thermal field at the interface was shown in Fig. 5.7 
for a succinonitrile/SiCp system. Note that if V is slightly higher than V cr . trough 
formation is readily observable (Fig. 5.7, b). On the contrary, if V » V cr , 
engulfment is very rapid, and no significant trough is observed experimentally 
(Fig. 5. 7, a). A comparison between calculated and experimentally evaluated data 
for this system are given in Table 5.2. Note that the calculated critical velocity for 
n = 1.5 was 4 mm/s, while the experimental data showed a critical velocity 
between 9.2 and 21 mm/s. The predicted value is thus reasonably close to the 
experimental one. 

Experimental data obtained by Bolling and Cisse l 5 - 10 ! on a water/tungsten 
system given in Table 5.2, show that the critical velocity for a particle having a 
radius of 5 mm is approximately 1 .8 mm/s for a gradient of 10 °C/cm. Calculation 
of the V -- d curves for a 5 mm radius tungsten particle in water are shown in Fig. 
5.10. The data used for calculation are given in Table 5.1, Appendix of this 
chapter. From Fig. 5.10 a critical velocity of 0.8 mm/s is predicted for n = 1 .5, and 
of 0.1 mm/s for n = 1.8. Given the incertitude in the evaluation of n and Act q , the 
calculated value for n = 1.5 is in reasonable agreement with the experimental 
results. 


Table 5.2 Experimental and calculated data on the critical velocity for various 
systems. 


System 

Particle 
dia., mm 

Experimental 

Calculated 

Interface 

Vcr 

mm/s 

Ref. 

Vcr, mm/s 
n = 1.5 

u>J 

SCN/SiC 

8 

planar 

9.2 

■ 

4 

water/tunqsten 

10 

planar 

1.8 

RBI 

0.8 

|! < 1 

AI-2wt.%Mq/SiC 

50-150 

cellular 

400 

ISI 

80 


















Fig. 5.10 Evaluation of critical velocity for a 5 mm radius tungsten particle 
in water for two values of n in Eq.(5.2). 

The temperature gradient will influence both the critical velocity and the 
critical distance, according to Eq. (5.18). However, its influence is rather limited. 
The level of gradient effect can be assessed from Eq. 5.18. It can be seen on Fig. 
5.5 that as the gradient increases from 90.6 to 108 °C/cm, for n = 1.5, the critical 
velocity increases from 3 to 4 mm/s, while the critical distance decreases from 
10.24 to 8.74 mm. These trends agree with the influence of gradient proposed in 
the numerical model of Potschke and Rogge I 5 - 6 ). They suggest that for a given 
particle radius the influence of the gradient is given by a relationship V cr = cons. • 
G 0 - 5 . This kind of dependency is in line with the experimental measurements 
reported by Korber et al. f 5 - 11 ! on the water -- latex particles system. 


It has been demonstrated by many authors [ 5 . 1 , 5 . 6 , 5 . 7 , 5 . 8 ] that the critical 
velocity decreases as the particle radius increases, following a power law. 
Calculations with the proposed model presented in Fig. 5.11 also confirm this 
type of dependency. 
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Fig. 5.1 1 . Relation between critical velocity and radius of particle 
for the succinonitrile/SiCp system, using n=1.5. 


B. Case 2: Kp/Km < 1 

According to the proponents of the "thermal conductivity criterion" model 
[5.4,5. 1 2 ] i W hen p = K y^ L < 1 . a bump rather than a depression should form at 

the liquid/solid interface. Indeed, a qualitative analysis of the thermal field by 
Bolling and Cisse t 5 - 10 !, as well as the analytical model proposed by Stefanescu 
et al. I 5 - 7 - 5 - 8 ], suggest the formation of such bumps on the interface. However, 
these analyses ignore the latent heat evolved during solidification. This heat can 
be "trapped" between the particle and the interface, because of the lower thermal 
conductivity of the particle. Thus the interface will be depressed, and a 
depression rather than a bump may form. However, when Kp/K|_ « 1, a bump 
can form. 

Some experimental data reported in the literature for systems that have m 
< 1 are given in Table 5.2. Unfortunately, the data for the two systems based on 
aluminum matrix are not for planar interface but for columnar or cellular interface. 
It is expected that the morphology of the interface will affect the critical velocity 
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[5-7] j although no quantitative treatment has been proposed so far. For the Al/SiC 
system the experiment [ 5 - 7 l showed that most of the particles were engulfed at V 
= 400 mm/s, but some were still pushed. Thus, this velocity can be considered 
slightly smaller than the critical velocity. The calculated value of 80 mm/s, as 
shown in Fig. 5.12, is reasonably close to the experiment. 
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Particle - interface distance, m 


Fig. 5.12. Influence of assumption made for n in Eq.(5.2) on the calculated 
critical distance for the aluminum/SiC system. Particle radius is 25 pm. 


5.5 Summary 

Through combined thermal field and force field calculation a model for 
prediction of particle behavior at the liquid -- solid interface was developed. The 
model is based on the hypothesis that both thermal equilibrium and force 
equilibrium must be achieved at the interface during steady state pushing of the 
particle when the particle is pushed at a distance higher than the atomic distance. 
If thermal equilibrium cannot be achieved, the particle can still be pushed, but the 
separation distance between the particle and the interface is of the order of the 
atomic distance. Two basic equations were derived. They were used for the 
calculation of the thermal and force equilibrium distance as a function of the 
solidification velocity. By coupling these equations a unique value can be found 
for the critical velocity. 
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The thermal field equation derived predicts the formation of bumps and 
troughs on the melt interface behind the particles. However, it demonstrates that 
the condition Kp/K|_ < 1 is not sufficient for bump formation. Limited effects of the 
thermal gradient on the critical velocity are also predicted. The conclusions of 
earlier researchers are thus confirmed. The occurrence of the bending of the 
interface was documented through experiments on transparent matrix materials 
(succinonitrile) containing SiC, and polystyrene particles. 

Model validation through calculation with data from the present 
experiments as well as from experiments performed by other researchers on 
water -- tungsten and aluminum -- SiC systems showed reasonable agreement. 
The model allows an approximate evaluation of the exponent used in the 
definition of the repulsive force based on surface energy. 

The present model is applicable to particulate metal matrix composites. It 
can be used to predict particle distribution in such composites solidifying with 
equiaxed structure, and inclusion distribution in castings. 

5.6 Appendix 2 


Table 1 Physical constants used in the calculation of the critical velocity. 


Physical 

constants 

' ' 

Succinonitrile - 
SiC 

Succinonitrile - 
Polystyrene 

Water - 
tungsten 

Aluminum - SiC 

Aluminum - 

ai 2 o 3 

a 0 (m) 

5 . 9 * 10- 10 

5 . 9 * 10" 10 

2 . 76 * 1 0 10 

1 . 43 * 1 0" 10 

1 . 43*1 O' 10 

h (kg m' 1 s" 1 ) 

0.0026 

0.0026 

0.001004 

0.0014 

0.0014 

K p (W m' 1 K' 1 ) 

85 

0.0794 

167.36 

85 

3.35 

K|_(W nr 1 K' 1 ) 

0.223 

0.223 

0.599 

93 

93 

r p (kg m‘ 3 ) 

3220 

1060 

19350 

3220 

3900 

r m (kg m' 3 ) 

986.7 

986.7 

998.2 

2357 

2357 

Cpp (J kg* 1 K' 1 ) 

1171.52 

- 

3.76 

- 

- 

C pm (J kg -1 K- 1 ) 

2026.95 

2026.95 

4183.0 

1079.5 

1079.5 

Ao 0 (N m" 1 ) * 

10 * 

10 * 

5 * 

8.27 171 

10 * 

Tsoak(K) 

353.0 

353.0 

298.0 

- 

- 

Tint (K) 

328.0 

328.0 

273.0 

- 

- 

G S (°C nr 1 ) 

8000 

8000 

8000 

3000 

- 

G L (“Cm" 1 ) 

10800 

10800 

10800 

9500 

- 

AH (J nr 3 ) 

4 . 6 * 1 0 7 

4 . 6 * 1 0 7 

3 . 33 * 1 0 8 

9 . 5 * 1 0 8 

9 . 5 * 1 0 8 


* Data are unavailable. These numbers were assumed to be reasonable. 




List of Symbols Chapter 5 


C pp 

Cpm 

Gl 

G S 

G*l 

Kp 

Kl 

Ks 

R 

T 

T m 

Tint 

Tp 

Tsoak 

V 

V C r 

V e 

3o 

d 

g 

m 

T| 

Pp 

Pm 

Ap 

Acr 0 

CT S p 

O s | 

CT|p 

AH 


Specific heat of the particle 
Specific heat of the matrix 

Temperature gradient in the liquid at distance L + R from the origin 

Temperature gradient in the solid at the solid/liquid interface 

Temperature gradient in the liquid film between the particle and the 

solid/liquid interface 

Thermal conductivity of particle 

Thermal conductivity of liquid matrix 

Thermal conductivity of solid matrix 

Particle radius 

Temperature 

Temperature in the matrix 

Temperature of the solid/liquid interface 

Temperature of the particle 

Soak temperature 

Interface velocity 

Critical velocity 

Equilibrium velocity 

Atomic distance 

Separation between the particle and the melt interface 
Acceleration due to gravity 
Thermal conductivity ratio (=Kp/K m ) 

Melt viscosity 
Particle density 
Matrix density 
Density difference 

Interfacial energy difference (= c sp - oi p - c s i ) 

Solid/particle interfacial energy 
Solid/liquid interfacial energy 
Liquid/particle interfacial energy 
Latent heat of fusion for the matrix 


Subscripts 
int: Interface 


I: Liquid m: Matrix p: Particle 
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